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Abstract 

We study the spectrum of the domain walls interpolating between dif- 
ferent chirally asymmetric vacua in supersymmetric QCD with the SU(3) 
gauge group and including 2 pairs of chiral matter multiplets in fundamental 
and anti- fundamental representations. For small enough masses m < m* ~ 
.286 Asqcd, there are two different domain wall solutions which are BPS- 
saturated and two types of "wallsome sphalerons". At m = m*, two BPS 
branches join together and, in the interval < m < m** 3.704 AsqcDj 
BPS equations have no solutions but there are solutions to the equations of 
motion describing a non-BPS domain wall and a sphaleron. For m > m**, 
there are no solutions whatsoever. 



1 Introduction 



Supersymmetric QCD is the theory involving a gauge vector supermultiplet V and 
some number of chiral matter supermultiplets. The models of this class attracted 
attention of theorists since the beginning of the eighties and many interesting and 
non-trivial results concerning their non-perturbative dynamics have been obtained 
The dynamics depends in an essential way on the gauge group, the matter 
content, the masses of the matter fields and their Yukawa couplings. 

The most simple in some sense variant of the model is based on the SU(N) 
gauge group and involves N — 1 pairs of chiral matter supermultiplets Si a , S" in 
the fundamental and anti-fundamental representations of the gauge group with a 
common mass m. The lagrangian of the model reads 



£ = f^Tr J d 2 6 W 2 + H.c.^ + \j d2ed2e S i eVS i 



+ - A J d 2 8d 2 6 s'^sl - y (/ d 2 6 S& + H.c. 



(1) 



color and Lorentz indices are suppressed. In this case, the gauge symmetry is broken 
completely and the theory involves a discrete set of vacuum states. The presence 
of N chirally asymmetric states has been known for a long time. They are best 
seen in the weak coupling limit m <C Asqcd where the chirally asymmetric states 
involve large vacuum expectation values of squark fields (sj) ^> Asqcd and the 
low energy dynamics of the model is described in terms of the colorless composite 
fields M.ij = ZS'iSj. The effective lagrangian presents a Wess-Zumino model with 
the superpotential 

^ = 4^-?^ (2) 
3 det.M 2 v ; 

The second term in Eq. (0) comes directly from the lagrangian (P and the first term 
is generated dynamically by instantons. Assuming Aiij = X 2 5ij and solving the 
equation dW/dx = (x is the scalar component of the superfield X), we find iV 
asymmetric vacua 

<X>*=(P^) (3) 
\3 m J 

(the vacua "k" and "k + N" have the same value of the moduli (x 2 )k and are 
physically equivalent). These vacua are characterized by a finite gluino condensate 

(Tr A 2 ) fe = 8n 2 m( X 2 ) k (4) 

It was noted recently that on top of (0) also a chirally symmetric vacuum with 
the zero value of the condensate exists. It cannot be detected in the framework of 
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Eq.(0) which was derived assuming that the scalar v.e.v. and the gluino condensate 
are nonzero and large, but is clearly seen if writing down the effective lagrangian 
due to Taylor, Veneziano, and Yankielowicz (TVY) || involving also the composite 
field 

$ 3 = -ArTr W 2 (5) 
32tt 2 v ; 

The corresponding superpotential reads Q 



2 n 

W = -$ 3 
3 



$ 3 detM 

^ n a 2JV+1 
IV SQCD 



~ yTr M (6) 



The presence of different degenerate physical vacua in the theory implies the 
existence of domain walls — static field configurations depending only on one spatial 
coordinate (z) which interpolate between one of the vacua at z = — oo and another 
one at z = oo and minimizing the energy functional. As was shown in in many 
cases the energy density of these walls can be found exactly due to the fact that the 
walls present the BPS-saturated states. 

The energy density of a BPS-saturated wall in SQCD with SU(N) gauge group 
satisfies a relation || 

A (Tr A 2 )oo - (Tr A 2 )^ (7) 



8tt 2 

where the subscript ±oo marks the values of the gluino condensate at spatial in- 
finities. 0. Bearing Eqs. (|7|, |j) in mind , the energy densities of the BPS walls 

are 

„ (im N - 1 V"' 

for the real walls and 

7T 

e c = 2e r sin — (9) 

for the complex walls. The RHS of Eqs.(|7]-[|) presents an absolute lower bound for 
the energy of any field configuration interpolating between different vacua. 

The relation (0) is valid assuming that the wall exists and is BPS-saturated. 
However, whether such a BPS-saturated domain wall exists or not is a non-trivial 
dynamic question which can be answered only in a specific study of a particular 
theory in interest. This question has already been studied in our previous works 
B B 0- m particular, in |5|, [7|, § the simplest model of the class (|l|) with 
N c = 2, Nf = 1 was analyzed. The results are the following: 

1 The factor 2/3 in Eq.(||) and the corresponding factor in Eq.(^) match the chosen (by historical 
reasons) normalization factor in the definition (pj) of 

2 A relation of this kind can be derived also for other variants of the theory involving exotic 
groups and more complicated matter content, but in general case the energy of a BPS wall and 
the gluino condensate are not related so directly H 
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1. For any value of the mass of the matter fields m, there are domain walls 
interpolating between a chirally asymmetric and the chirally symmetric vacua 
(we call them real walls). They are BPS - saturated. 

2. There are also complex BPS solutions interpolating between different chirally 
asymmetric vacua. But they exist only if the mass is small enough m < m* = 
4.67059 . . .Asqcd- When m > m*, BPS walls are absent. 

3. In a narrow range of masses < m < m** « 4.83 Asqcd, complex domain 
walls still exist, but they are not BPS saturated anymore. At m > m**, there 
are no such walls whatsoever. 

In Ref.0, we studied the problem of existence of BPS-saturated domain walls in 
the model ([!]) with N > 3. The results are basically the same as for N = 2: the 
real walls exist for any m and are BPS-saturated, and there are two complex BPS 
branches which exist in a limited range m < m*. The value of m* goes down with 
N: m* = .28604 . . . Asqcd for SU(3), m* = .07539 ... A SQC d for 577(4), etc. For 
large N, m*(N) oc N' 3 . 

The results concerning the BPS walls were obtained by solving numerically the 
first order BPS equations 

d z( j) = e iS dW/d0, d zX = e l5 dW/dx (10) 

associated with the TVY lagrangian. The phase S depends on particular vacua 
between which the wall interpolates (see Refs.|| [10], || for details). 

To study the spectrum of the domain walls which are not BPS-saturated, one 
has to solve the equations of motion which are of the second order, and, technically, 
the problem is a little bit more involved. We did it earlier for iV = 2 ||. This paper 
is devoted to the numerical solution of the equations of motion for the 577(3) gauge 
group. 



2 Solving equations of motion 

The scalar potential corresponding to the superpotential 



is 





2 


dW 


2 




+ 


dx 


= 4 







2 ln{0V (Ar_1) } +(N-lf 



mx 



3X 



'ID 



(from now on we set Asqcd = 1)- The potential (|TT|) has N + 1 degenerate minima. 
One of them is chirally symmetric: = x — 0. There are also N chirally asymmetric 
vacua with (x)k given in Eq.(Q) and 



2i(JV-l)-n-fc 

e sjv 



(12) 
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To study the domain wall configurations, we should add to the potential ( [TT]) 
the kinetic term which we choose in the simplest possible form 



a™ = m 2 + \d X ? 

and solve the equations of motion with boundary conditions 

(0) o = R., 0(oo) = R^e- 2 ^ N -^ N 



(13) 



-oo 



x(-oo) = (x)o = P*, x(oo)=P*e 



ni/N 



(14) 



Thereby we are studying the walls interpolating between "adjacent" complex vacua. 
Actual calculations will be performed for N = 3 where all the vacua are adjacent. 
In principle, one could also study numerically the walls interpolating between the 
vacua k — and k = 2 for, say, SU (5) gauge group, etc. We expect the physical 
results for all such walls to be qualitatively the same. 

It is convenient to introduce the polar variables x = pe ta , <P = R e% ^ after which 
the equations of motion acquire the form 



R" - Rf3' 2 
R0" + 2R'p' 



p" - pa 12 



16R 5 
~3p^~ 



8R 3 [L(L + 3/2)+(3 2 + ] + (N - l) 2 
12R 3 P + + A(N - ifmR 2 sin(/3_) 

AT - 1 L + (N- 1 2 m p — - 

p V V 



AmR 2 cos(/3_ 



pa" + 2p'a 



(N 



(N-l) 



,8mR 3 



sin(/3_ 



(15) 



where L = 
conditions 



p 3p 
ln[R 3 p 2 ( N -V], (3 + = 3(3 + 2(N - l)a, /3_ = 3/3 — 2a, with the boundary 



p(— 00) = p(oo) = p*; R(—oo) = R(oo) = R*; 

2(N-l)ir 



a(-oo) = (3(-oo) = 0; a(oo) = tt/N; j3(oo) 



3N 



(16) 



When N = 2, the system fllBD is reduced to that studied in Ref.||. The system ( |15| 
involves one integral of motion 



T-U = R' 2 + p r2 + tt t p" + f?a r2 - 



,/2 



2,3/2 



,2/2 



4R\L 2 + p 2 + ) -{N-l) 



2 2 1 
m p + 



16R 6 8mR 3 



9p 2 



cos(/3_ 



const 



(17) 



In our case, const = due to boundary conditions fll6"l) . The phase space of the 
system ([15]) is 8-dimensional and a general Cauchy problem involves 8 initial condi- 
tions. The problem is simplified, however, when noting that the wall solution should 
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Figure 1: The ratio r) = R(Q)/Rq(0) for the solutions of the equations of motion as 
a function of mass for the SU(3) theory. The solid lines describe the BPS solutions 
and the dashed lines describe the non-BPS wall and the sphalerons. 



be symmetric with respect to its center. Let us seek for the solution centered at 
z = so that 



p(z)=p(-z), R(z) = R(-z), 
a(z) = tt/N - a{-z), (3{z) = -2(N - l)n/(3N) - (3{-z) (18) 

Indeed, one can be easily convinced that the Ansatz fll8[) goes through the equations 
(|15|). It is convenient to solve the equations ([15]) numerically on the half-interval 
from z = to z = oo. The symmetry (]T8|) dictates p'(0) = R'(0) = 0, a(0) = 
tt/(2N), (3(0) = — (N — l)n/(3N) which fixes 4 initial conditions. Four others 
satisfy the relation dl7|) . Thus, we are left with 3 free parameters, say, p(0), R(0), 
and f3'(0), which should be fitted so that the solution approach the complex vacuum 
in Eq.(f|) at z -> oo. 

All the solutions obtained with such a procedure are presented in Fig. [I] where 
the parameter R(0), one of the fitted initial conditions, is plotted as a function of 
mass (we normalized R(0) at its value Rq(0) = .918 . . . i?* for the upper BPS branch 
in the limit m — > 0; see Ref. for details). For small masses, there are several 
solutions. We obtain first of all the solutions studied in Ref. @ and describing the 
BPS-saturated domain walls, (solid lines in Fig. |T]). We find also two new solution 
branches drawn with the dashed lines in Fig. |I[ We see that, similarly to the BPS 
branches, two new dashed branches fuse together at some m = m** w 3.704 . No 
solution for the system ([15]) exist at m > m^. 

The picture is rather analogous to what we had for iV = 2 || and the physical 
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R(0) 
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R(0) 
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b) m* < m < m* 



c)m = m 



Figure 2: Illustrative profiles of the energy functional vs. R(0) 



interpretation is similar. Let us assume first that m < and draw the energy 
functional E for field configurations with wall boundary conditions minimized over 
all parameters except the value of R(0) which is kept fixed (see Fig. ^|a). For very 
small -R(O), our configuration nearly passes the chirally symmetric minimum and 
the minimum of the energy corresponds to two widely separated real walls. Thus 
E[R(0) = 0] = 2e r with e r given in Eq.(^). Two minima in Fig. 0a correspond to 
BPS solutions with the energy e c = y3e r . They are separated by an energy barrier. 
The top of this barrier (actually, this is a saddle point with only one unstable mode 
corresponding to R(0), in other words — a sphaleron) is a solution described by the 
upper dashed line in Fig.|I|. The lower dashed line corresponds to the local maximum 
on the energy barrier separating the lower BPS branch and the configuration of two 
distant real walls at R(0) = 0. f] 

At m — m*, two BPS minima fuse together and the energy barrier separating 
them disappears. The upper sphaleron branch coincides with the BPS solution at 
this point. When m is increased above m#, the former BPS minimum is still a 
minimum of the energy functional, but its energy is now above the BPS bound (see 
Fig.|2|b). The corresponding solution is described by the analytic continuation of the 
upper sphaleron branch. The lower dashed branch in the region m* < m < m** is 

3 In contrast to the case N = 2, the existence of such a barrier could not be established from 
the BPS spectrum alone. The matter is that, while 2e r = e c for N = 2 and the presence of the 
maximum is guaranteed by the Roll theorem, 2e r > e c in our case and one could in principle 
imagine a situation where E falls down monotonically when R(0) is increased from zero up to its 
value at the lower BPS branch. As will be discussed later, our numerical results are not good 
enough to make a rigid statement on the form of the function E[R(0)] at small R(0) in the region 
to < to*, but they strongly suggest that the energy barrier ( though a tiny one) is present. 
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Figure 3: The ratio n(m) = E/e c — 1 for the non-BPS wall and for the lower 
sphaleron as a function of mass. 

still a sphaleron. At the second critical point m = m**, the picture is changed again 
(see Fig.|2]c). The local maximum and the local minimum fuse together and the only 
one remaining stationary point does not correspond to an extremum of the energy 
functional anymore. At larger masses, no non-trivial stationary points are left. 

Our findings are illustrated in Figs. |3], ^ where the energies of the non-BPS 
wall and of the lower sphaleron branch are plotted as a function of m. In Fig. |3], 
the ratios of the energies of both branches to the BPS bound (H) are plotted. The 
lower line in Fig. |3| corresponds to the stable wall solution and the upper line to 
the sphaleron branch. For almost all m < m**, the wall solution is globally stable. 
When m** — m becomes very small, it is stable only locally: we see that, at m = m** 
where two branches are fused together, their energy exceeds slightly the energy of 
two real walls 2e r . 

In Fig. ^ the sphaleron energy is redrawn in logarithmic scale in the units of 
2e r . Unfortunately, we do not have good numerical data at m < .7 — our relative 
uncertainty becomes large. We see, however, that the logarithmic plot in Fig. |] is 
pretty much linear, the best fit is 

«'(m) = 5.49 • 10" 7 exp{2.11m} (19) 

The fit (|T9|) cannot be valid for very small masses: we expect «/(0) = which means 
that the straight line in Fig. § should bend down at small enough masses due to a 
preexponential factor ~ m a which we cannot determine from our data. Anyway, it 
is seen from Fig. f| that In k' does not "want" to hit infinity at a finite mass, and 
we assume that it does not ( though we cannot exclude it as a logical possibility). 
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Figure 4: Logarithmic plot for the ratio k'{w) = Ei OWCT sp haieron/(2e r ) — 1 vs. mass. 

In terms of Fig. |2], that means that the energy barrier on the left (for illustrative 
purposes, it is very much exageratted ) is present for all masses. 

Finally, we present for illustration the profile R(z) for the lower sphaleron branch 
at m = .7. As was expected, it resembles very much a combination of two separate 
real domain walls. 

3 Discussion. 

Our main conclusion is that, besides the critical mass beyond which BPS solu- 
tions disappear, also a second critical mass m** exists beyond which no complex wall 
solution can be found whatsoever. This was the case for SU{2) || and, as we see 
now, this is also the case for SU(3). Seemingly, the same situation holds for any N. 
That means in particular that no domain walls connecting different chiraly asym- 
metric vacua are left in the pure supersymmetric Yang-Mills theory corresponding 
to the limit m — > oo, and only the real domain walls connecting the chirally symmet- 
ric and a chirally asymmetric vacuum states survive in this limit. That contradicts 
an assumption of Ref.[T]J that it is complex rather than real domain walls which 
are present in the pure SYM theory (Witten discussed them in the context of brane 
dynamics) . 

One has to make a reservation here: our result was obtained in the framework of 
the TVY effective lagrangian (||) whose status [in contrast to that of the lagrangian 
(0)] is not absolutely clear: the field $ describes heavy degrees of freedom (viz. a 
scalar glueball and its superparnter) which are not nicely separated from all the rest. 
However, the form of the superpotential (P) and hence the form of the lagrangian 
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Figure 5: Profile R(z) for the lower sphaleron branch at m — .7 



for static field configurations is rigidly dictated by symmetry considerations; the 
uncertainty involves only kinetic terms. It is reasonable to assume that, as far as 
the vacuum structure of the theory is concerned (but not e.g. the excitation spec- 
trum — see Ref. || for detailed discussion), the effective TVY potential @ can be 
trusted. A recent argument against using Eq. (§) that the chirally symmetric phase 
whose existence follows from the TVY lagrangian does not fulfill certain discrete 



anomaly matching conditions |12| is probably not sufficient. First, it assumes that 
the excitation spectrum in the symmetric phase is the same as it appears in the TVY 
lagrangian which is not justified. Second, it was argued recently that the TVY la- 
grangian describes actually all the relevant symmetries of the underlying theory and 
the absence of the anomaly matching is in a sense an "optical illusion" |13| . 

The main distinction of the SU(3) case considered here compared to the SU(2) 
theory is that the values of two critical values are rather different (in SU(2) case 
they were pretty close: m* = 4.67059 . . . and m** « 4.83). This is due to the fact 
that the energy of the complex BPS wall e c is less in this case than the energy of 
two real walls 2e r . When we increase the mass and go above m*, the energy of the 
wall first has to rise from e c to 2e r . Only then the complex domain wall "bound 
state" can break apart into its "constituents" , the real walls. Q 

One can expect that m, and m** differ more and more as N grows. A tentative 
guess is that m** is roughly iV-independent (to be compared with m*(N) oc iV~ 3 ). 
Of course, that can be confirmed or disproved by only actual numerical study. Note, 
however, that numerical calculations become more and more difficult as N grows 



4 Actually, as we have seen, the energy barrier in Fig. || does not allow the complex wall to 
break apart until its energy goes a little bit above the limit 2e r . 
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- the instabilities characterized by eigenvalues of the Jacobi matrix of the system 
(|15|) near the minima 0) grow as N 2 . 
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